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LONG-TERM  GOALS 

Wind-generated  waves  play  a  prominent  role  at  the  interfaces  of  the  ocean  with  the  atmosphere,  land 
and  solid  Earth.  Waves  also  define  in  many  ways  the  appearance  of  the  ocean  seen  by  remote-sensing 
instruments.  Beyond  these  geophysical  aspects,  waves  also  affect  human  activities  at  sea  and  on  the 
coast.  The  long-term  goals  of  this  research  are  to  obtain  a  better  understanding  of  the  physical  processes 
that  affect  ocean  surface  waves  and  their  interactions  with  ocean  currents  and  turbulence,  the 
atmosphere,  seismic  waves,  sediments  and  remote  sensing  systems,  and  to  improve  our  forecasting  and 
hindcasting  capacity  of  these  phenomena  from  the  global  ocean  to  the  nearshore  scale. 

OBJECTIVES 

•  Observe  and  parameterize  the  dissipation  of  ocean  waves  due  to  breaking,  wind-wave 
interactions,  or  bottom  friction 

•  Advance  spectral  wave  modeling  at  all  (global  to  beach)  scales  in  a  unified  framework,  in  terms 
of  parameterization  and  numerical  developments 

•  Help  the  application  of  wave  models  to  new  problems  (upper  ocean  mixing  and  surface  drift,  use 
of  seismic  noise  data,  air-sea  gas  exchange  ...)  and  use  these  applications  for  feedback  on  the 
wave  model  quality 

APPROACH 

By  combining  theoretical  advances  with  numerical  models,  remote  sensing  and  field  observations,  we 
investigate  the  physical  processes  that  affect  wind-generated  ocean  gravity  waves.  The  various 
dissipative  processes  that  contribute  to  the  spectral  wave  evolution  are  isolated  by  considering 
geophysical  situations  in  which  they  are  dominant:  the  long-distance  swell  propagation  in  the  case  of 
air-sea  friction,  the  evolution  of  swells  on  shallow  continent  shelves  in  the  case  of  bottom  friction,  the 
energy  level  in  the  spectral  tail  in  the  case  of  cumulative  breaking  effects,  and  the  breaking  statistics  of 
waves.  These  require  the  acquisition  of  new  data  using  stereo- video  techniques,  for  the  spectral  levels 
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of  waves  of  1  to  10  m  wavelength,  and  the  statistics  of  whitecaps.  The  full  model  is  then  confronted  to  a 
wide  range  of  observations  starting  from  global  altimeter,  SAR  and  buoy  data.  Alvise  Benetazzo 
(ISMAR  Venezia)  is  performing  the  calibration  of  the  stereo  system  and  the  reconstruction  of  sea 
surface  geometries,  that  expertise  has  been  transferred  to  Ifremer  with  graduate  student  Fabien  Leckler 
and  developments  at  the  two  institutions  are  now  allowing  reliable  analysis  of  data  acquired  in  201 1. 
The  spectral  analysis  and  whitecap  detection  is  performed  at  Ifremer,  under  the  supervision  of  Fabrice 
Ardhuin.  All  the  wave  modeling  effort  at  Ifremer  (theory,  parameterization  and  calibration)  is 
performed  by  Fabrice  Ardhuin  and  Fabien  Leckler. 

WORK  COMPLETED 

Wave  model  parameterizations  for  wind  wave  evolution  The  main  challenge  faced  for  the  energy 
balance  of  wind  seas  is  the  non-local  nature  of  all  terms  in  the  spectral  domain.  This  is  well  known  for 
the  non-linear  4-wave  interaction  term  (Hasselmann,  1962),  but  it  is  becoming  also  increasingly  clear 
for  the  wind  input  and  the  dissipation  terms.  For  the  wind  input,  the  well-known  quasi-linear  effect 
(Fabrikant,  1976;  Janssen,  1991)  is  a  positive  feedback  of  the  short  waves  in  the  spectrum  on  the  growth 
of  all  spectral  components,  whereas  the  sheltering  effect  discussed  by  Kara  and  Belcher  (2004),  is  a 
negative  feedback  term  that  affects  the  short  wave  growth  based  on  the  long  wave  growth.  The 
combination  of  the  two  effects  leads  to  behaviors  that  are  not  fully  understood,  with  clear  impacts  on 
the  variability  of  wind  stress  and  mean  square  slopes. 

Similarly,  there  is  growing  evendence  for  the  generation  of  both  longer  and  shorter  wave  components 
from  the  breaking  of  waves  of  intermediate  scale  (Yurovskaya  et  ah,  2008;  Tulin  and  WAseda,  1999; 
Caudal,  2012),  in  addition  to  the  “cumulative”  effect  of  short  wave  dissipation  due  to  large  scale 
breakers,  already  introduced  in  source  term  parameterizations  (Banner  and  Morison,  2010;  Ardhuin 
et  ah,  2010;  Rogers  et  ah,  2010).  On  top  of  that,  it  is  also  expected  that  the  breaking  of  waves  changes 
the  wind  profile  and  wind-wave  growth  term  (Reul  et  ah,  2008),  an  effect  introduced  in  the 
parameterization  by  Banner  and  Morison  (2010). 

Finally,  it  is  well  known  that  the  DIA  parameterization  for  non-linear  interactions  leads  to  biases  in 
spectral  shape  and  integrated  parameters  (Banner  and  Young,  1994;  Ardhuin  et  ah,  2007).  Using  a  very 
preliminary  test  of  the  multiple-DIA  from  (?),  Rollet  (2013)  showed  how  important  the  directional 
spreading  was  for  far-held  swells  from  a  storm.  Indeed  the  settings  of  the  MDIA  tuned  to  the  Tolman 
and  Chalikov  parameterization  produced  pathologically  narrow  directional  spectra  when  applied  with 
the  TEST451  parameterization  from  Rascle  and  Ardhuin  (2013).  This  was  a  good  candidate  to  reveal 
that  the  shape  of  the  swell  helds  in  extra  tropical  storms  is  strongly  inhuenced  by  the  directional 
spreading  of  the  spectra  within  the  storm,  with  different  results  given  by  different  parameterizations. 
This  is  illustrated  in  hgure  1.  Future  work  with  MDIA  will  require  a  tuning  of  the  interaction 
coefficients. 

There  are  thus  many  effects  that  need  to  be  considered,  with  strong  non-linear  interactions  between 
them,  and  it  is  not  clear  how  to  transform  each  type  of  measurement  into  a  clear  constraint  for  the 
functional  form  and  magnitude  of  these  effects.  We  have  begun  to  test  several  aspects  and  there  is  not 
clear  positive  trade  off  when  adding  more  effects  into  the  TEST451  parameterization  (?).  This  work  has 
been  partly  delayed  due  to  the  difficulty  of  hiring  personnel  at  Ifremer  in  EY2013.  This  line  of 
investigation  will  be  continued. 
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(a)  derived  from  Envisat  ASAR  (b)  Model  with  +DIA  (c)  Model  with  MDIA 


Figure  1:  Observed  and  modelled  maps  of  significant  swell  height,  one  day  after  propagation  of  swell 
from  April  3,  2009,  in  the  North  Pacific.  The  MDIA  parameterization  of  4-wave  interactions  is  used 
with  the  settings  adjusted  to  the  Tolman  and  Chalikov  wind-wave  growth  and  dissipation  and  should 
be  retuned  for  the  TEST451  parameterization  used  here.  Picture  adapted  from  Rollet  (2013). 


Bottom  friction  and  coastal  model  numerics  Bottom  friction:  The  parameterization  by  Ardhuin  et  al. 
(2003)  based  on  data  of  the  1999  SHOaling  Waves  Experiment  (SHOWEX)  has  been  thoroughly  tested 
for  the  North  Carolina  shelf  and  Erench  Atlantic  coast.  Adjustements  included  the  setting  of  maximum 
roughness  values  for  rock  bottoms. 

Hindcasts  and  wind  speed  adjustements  The  time- varying  biases  in  ECMWE  operational  analyses 
and  CESR  hindcasts  require  some  specific  adjustments  of  model  parameters.  Eor  recent  years,  the 
higher  level  of  wind  speed  percentiles  above  80  require  a  specific  reduction  in  wind  speed  for  the  years 
2007  to  2012.  We  have  thus  corrected  oour  hindcasts  for  these  years  with  a  reduced  BETAMAX 
parameter  from  1.33  to  1.15.  Similarly,  the  value  used  with  ECMWE  operational  anlyses  was  reduced 
from  1.52  for  2009  to  1.45  for  the  year  2011,  and  further  reductions  are  expected  with  the  increase  in 
resolution  of  the  model  and  change  of  WAM  model  parameterization. 

Data  analysis 

Acoustic  noise:  Ardhuin  et  al.  (2013)  have  used  the  variability  of  ocean  bottom  noise  records  at 
frequencies  below  1  Hz  to  validate  the  changing  directional  properties  of  ocean  waves  as  modeled  with 
different  source  term  parameterizations  parameterizations.  Stereo-video  data:  The  quality  of  3D 
reconstructions  is  now  well  established  as  discussed  in  Eedele  et  al.  (2013)  based  on  data  acquired  in 
2009.  Recordings  from  the  201 1  experiment  in  Katsiveli  (Black  Sea),  have  finally  been  processed  to  a 
satisfactory  degree  of  quality,  thanks  ot  to  the  contribution  from  Alvise  Benetazzo  who  developped  a 
self-calibration  procedure  that  allowed  to  correct  the  viewing  angles  of  cameras. 

Implementation  in  WAVEWATCH  III 

Several  items  have  made  their  way  to  the  trunk  of  the  NCEP  subversion  server  over  EY2013 

•  an  upgrade  of  the  multi-grid  system  allowing  two-way  nesting  between  regular  grids  and 
triangle-based  meshes. 

•  an  interface  to  the  lOS  tidal  analysis  software,  allowing  water  levels  and  currents  to  be  defined 
from  tidal  constituents. 

•  estimation  of  whitecap  coverage  and  foam  thickness  from  the  dissipation  parameterization  as 
described  in  Eeckler  et  al.  (2013) 
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While  others  are  still  being  tested: 

•  a  data  assimilation  procedure  based  on  SAR-derived  swell  fields. 

•  The  Eldeberky  and  Battjes  triad  parameterization 

•  a  parametric  source  of  free  infragravity  wave  energy  at  the  shorelines. 

•  an  evolution  of  wind  wave  generation  and  dissipation,  taking  into  account  modulation  effects  and 
short  wave  generation  by  breaking  waves. 

RESULTS 

Most  of  the  development  work  this  year  has  focused  on  coastal  applications,  in  particular  with 
numerical  aspects  associated  to  unstructured  grid,  for  which  the  WAVEWATCH  III  implementation  has 
proven  very  robust,  and  calibration  work  on  bottom  friction.  This  follows  previous  works  on  the  effects 
of  currents,  described  in  ?.  Bottom  friction.  Using  the  adjusted  SHOWEX  bottom  friction  (Ardhuin 
et  ah,  2003),  a  20-year  hindcast  of  coastal  sea  states  around  Erance,  using  tidal  currents  and  elevations, 
and  a  map  of  sediment  types.  Roland  and  Ardhuin  (tted)  have  shown  that  this  is  the  most  accurate 
hindcast  produced  so  far  for  that  region,  and  the  variability  of  bottom  types  does  correspond  to  different 
wave  attenuation  patterns  recorded  by  buoys.  This  is  illustrated  in  figure  2. 

Underwater  acoustics  Another  area  of  activity  was  the  application  of  numerical  wave  models  to  the 
analysis  of  underwater  sound  at  frequencies  below  1  Hz.  Ardhuin  et  al.  (2013)  showed  that  the  latest 
parameterization  by  Rascle  and  Ardhuin  (2013)  was  particularly  well  suited  to  reproduce  the  variability 
of  the  acoustic  noise  level  induced  by  the  directional  wave  distribution,  at  acoustic  frequencies  up  to 
0.8  Hz  in  the  case  of  data  recorded  north  of  Hawaii.  At  higher  frequencies,  the  model  misses  the 
dynamics  of  the  recorded  signal,  pointing  to  the  missing  impact  of  wave  breaking  on  the  wave 
spectrum,  as  already  suggested  by  Duennebier  et  al.  (2012),  and  consistent  with  radar  backscatter 
studies  (Yurovskaya  et  al.,  2008). 

IMPACT  AND  APPLICATIONS 

The  combined  use  of  underwater  acoustic  data  and  numerical  wave  models  offers  new  opportunities  for 
investigating  the  directional  wave  spectrum,  up  to  very  high  frequencies.  We  are  aware  of  the  work 
done  in  San  Diego  by  W.  Earrell  and  W.  Munk  on  this  topic  and  have  started  exchanging  ideas  about  it. 
W.  Earrell’s  early  comments  on  our  acoustic  noise  modeling  paper  by  Ardhuin  et  al.  (2013)  was 
essential  for  resolving  discrepancies  in  previous  analyses.  We  are  currently  processing  the  2011 
Katsiveli  data  to  reveal  the  shape  of  1  to  5  m  wavelength  waves,  without  directional  ambiguity,  in 
connection  with  that  particular  acoustic  work. 

National  Security  Improving  wave  forecasts  are  relevant  to  a  variety  of  denfence  applications.  The 
most  dramatic  improvement  brought  in  the  operational  models  by  the  preasent  work  is  an  improved 
representation  of  swells  which  are  most  relevant  for  amphibious  operations. 

Quality  of  Life  The  transport  of  contaminants  in  the  nearshore  ocean  is  largely  driven  by  waves.  The 
capability  and  understanding  of  this  driving  process  in  three  dimensions  will  certainly  lead  to  improved 
water  quality  models. 


4 


52 
51 
50 
49 
48 
47 
46 
45 
44 
43 

-6-4-2  0  2  4 


.^^62059  (b) 

-^62069 

\ 

? 

!62067 

.  1  mm 

.  0.1mm 

\  0.01mm 


I 


Figure  2:  (a)  Map  of  sediment  median  diameter  and  (b)  mean  difference  in  significant  wave  height 
(in  meters)  over  the  month  of  February  2010  between  a  model  run  using  the  ’JON-  SWAP’  bottom 
friction  parameterization  and  another  using  the  ’SHOWEX’  parameterization  with  a  constant 
Nikuradse  roughness  length  of  12  cm  for  rocks.  Inset  is  a  zoom  on  the  region  around  Yeu  and 
Noirmoutier  islands  where  the  impact  of  this  friction  is  very  clear,  (c)-(e),  Time  series  of  observed 
and  modeled  significant  wave  height  at  several  buoys  using  the  JON SWAP  (blue  diamonds)  or 
SHOWEX  (red  triangles)  parameterizations  for  bottom  friction,  compared  to  hourly  buoy 
measurements  (solid  line).  The  location  of  these  buoys  is  indicated  by  stars  in  (b).  Figure  adapted 

from  Roland  and  Ardhuin  (tied). 


TRANSITIONS 

As  mentionned  above,  operational  wave  models  at  NOAA/NCEP  have  been  switehed  from  the 
paramterization  by  Tolman  and  Chalikov  (1996)  to  the  TEST451,  only  3  months  after  the 
parameterization  had  been  adjusted.  This  shows  the  amazing  eapability  brought  about  by  the  new  way 
of  developing  model  frameworks  at  NCEP  using  outside  eontributions  on  the  same  subversion  server. 
This  happened  during  PY2012.  Sinee  that  time,  many  adjustments  and  bug  eorreetions  have  been 
disseminated,  using  the  same  tools,  and  have  made  their  way  into  the  version  4.12  of  the 
WAVEWATCH  III  eode,  now  publieally  available  on  request  from  NCEP  We  will  eontinue  helping 
NOAA/NCEP  and  others  in  testing  and  implementing  these  parameterizations.  We  are  also  eomitted  to 
support  outside  users,  with  a  training  eourse  in  Brest,  Eranee,  whieh  is  open  to  all  interested  parties,  in 
eoordination  with  the  University  of  Maryland  and  NCEP  efforts.  The  first  training  will  take  plaee  on  the 
week  of  November  4,  2013  and  partieipants  from  all  over  the  world  will  attend,  ineluding  people  from 
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Deltares  (NL),  and  U.S.  universities.  A  seeond  edition  is  seheduled  for  Mareh  2014. 

RELATED  PROJECTS 

The  present  “Oeean  Waves  Dissipation  and  speetral  Balanee”  (WAVE-DB)  shares  many  of  the 
objeetives  of  the  the  Integrated  Oeean  Waves  for  Geophysieal  and  other  Applieations  (lOWAGA) 
projeet,  funded  by  the  European  Researeh  Couneil.  As  a  result,  results  from  both  projeets  are  reported 
on  the  same  web  pages,  where  the  eontribution  from  eaeh  is  elearly  identified.  Whereas  WAVE-DB  is 
foeused  on  the  development  of  stereo-video  teehniques  and  numerieal  wave  modeling,  lOWAGA  allows 
allows  a  broader  perspeetive  with  work  on  remote  sensing  and  seismie  noise,  whieh  allow  a  more 
informed  ealibration  of  the  numerieal  wave  model.  Einally,  the  WAVE-DB  aetivity  is  also  benefiting 
from  the  GEOBWAVE  projeet,  funded  by  the  European  Spaee  Ageney  and  the  Ereneh  Spaee  Ageney  to 
faeilitate  the  use  of  satellite  remote  sensing  data  of  oeean  waves. 
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